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In the Relativistic Mean Field Approximation three quarks in baryons from the lowest octet and 
the decuplet are bound by the self-consistent chiral field, and there are additional quark-antiquark 
pairs whose wave function also follows from the mean field. We present a generating functional for 
the 3-quark, 5-quark, 7-quark ... wave functions inside the octet, decuplet and antidecuplet baryons 
treated in a universal and compact way. The 3-quark components have the 517 (6)-symmetric wave 
functions but with specific relativistic corrections which are generally not small. In particular, the 
normalization of the 5-quark component in the nucleon is about 50% of the 3-quark component. 
We give explicitly the 5-quark wave functions of the nucleon and of the exotic + . We develop a 
1^-^ , formalism how to compute observables related to the 3- and 5-quark Fock components of baryons, 

("*^ ■ and apply it to estimate the + width which turns out to be very small, 2-4 MeV, although with a 

f"^ ' large uncertainty. 
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PACS numbers: 12.38.-t, 12.39.-x, 12.39.Dc, 14.20-c 



I. INTRODUCTION 



£S| ' Were the chiral symmetry of the QCD lagrangian not broken spontaneously, the nucleon would be either nearly 
massless or degenerate with its chiral partner, iV(1535, i ). Both alternatives are many hundreds of MeV away from 
reality, which serves as one of the most spectacular indications that chiral symmetry is spontaneously broken. It also 
serves as a warning that if we disregard the effects of the spontaneous chiral symmetry breaking we shall get nowhere 

■ in understanding light baryons. 

Spontaneous chiral symmetry breaking implies that at the microscopic level of QCD nearly massless u, d, s quarks 

■ gain a dynamical momentum-dependent mass M(p) with M Ut d(0) ~ 350 MeV. A probable mechanism [1] of how it 
happens is provided by instantons - large fluctuations of the gluon field in the vacuum. The resulting massive quarks 
are usually called the constituent quarks; they necessarily, as a consequence of chiral symmetry, have to interact 

O ' with the (pseudo) Goldstone pion field, and actually very strongly: the dimensionless coupling constant is about 
A/(0)/F 7r w 4. The corresponding low-energy interaction lagrangian is written below, in Section II. It implies that 

J^" 1, inside baryons there is a strong chiral field. Generally speaking, the chiral field experiences quantum fluctuations; 

(-h | however, one may ask if it is reasonable to introduce the notion of a mean chiral field inside baryons. 

The mean field approach to bound states is usually justified by the large number of participants. The Thomas-Fermi 
. £h ' approximation to atoms is justified at large Z, and the shell model for nuclei is justified at large A. In baryons, the 
appropriate large parameter justifying the mean field approach would be the number of colors N c [2]. The number of 
colors being iV c = 3 in the real world, one may wonder how accurate is the mean-field picture. Theoretically speaking, 
there are two kind of corrections in 1/N C to the mean field. One kind is due to the high-frequency fluctuations of 
the chiral field about its mean-field value in a baryon. These are loop corrections and are additionally suppressed by 
factors of l/(27r). With the present precision, such corrections, typically of the order of 10%, can be ignored. The 
second type can be called kinematical: they are due to the rotations of the baryon mean field in ordinary and flavor 
spaces, and are not suppressed additionally. Such corrections are not small at N c = 3 (although they tend to zero in 
the academic limit N c — > oo) and should be taken into account exactly, if possible. 

In this paper, we adopt the view [3] that there is a self-consistent mean chiral field in baryons, which binds three 
massive constituent quarks, see Fig. 1. The binding appears to be rather tight; bound-state quarks are relativistic and 
their wave function has both the upper s-wave Dirac component and the lower p-wave Dirac component, see Section 
III. Simultaneously, the negative-energy Dirac sea of constituent quarks is distorted by the same mean field, leading to 
the presence of an indefinite number of additional quark-antiquark (QQ) pairs in baryons, see Fig. 2. Ordinary baryons 
are superpositions of 3Q, 5Q, 7Q... Fock components. This picture which we shall call the Relativistic Mean Field 
Approximation to baryons (or else the Chiral Quark Soliton Model where the word "soliton " is an alias of the mean 
field), leads, without any fitting parameters, to a reasonable quantitative description of the baryons properties [3, 4], 
including nucleon parton distributions at a low normalization point [5] and other baryon characteristics [6] . It should 
be stressed that the approximation supports full relativistic invariance and all symmetries following from QCD. 

We shall see that the normalization of the 5Q component in the nucleon is not small as compared to its 3Q 
component. The three-quark picture of a nucleon is an out-fashioned cartoon. It might do in popular lectures but 
professionals should explain why the spin carried by three quarks is three times less, and the nucleon tr-term is 
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four times bigger than in the naive 3Q picture [7]. Taking into account the QQ pairs in the nucleon explains these 
paradoxes [8, 9]. 
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FIG. 1: A schematic view of baryons in the Relativistic Mean 
Field Approximation. There are three "valence" quarks at a dis- 
crete energy level created by the mean field, and the negative- 
energy Dirac continuum distorted by the mean field, as com- 
pared to the free one. 



FIG. 2: Equivalent view of baryons in the same approxi- 
mation, where the distorted Dirac sea is presented as QQ 
pairs. The average number of QQ pairs is proportional 
to the amplitude squared of the mean field, times N c . 



The correct lowest baryons' quantum numbers arise as a result of the quantization of the rotation of the mean 
chiral field in the ordinary and in the flavor spaces [2]. If the mean chiral field is presented as a unitary 3x3 matrix 
V(x), the rotated field is 

U(x) = RV(x)rf (1) 

where R is an SU (3) rotation matrix; it can be parameterized by eight "Euler angles" as it is done, for example, in 
Appendix A. We shall for simplicity set the strange quark mass m s = 0; in this limit any rotated mean field (1) is, 
classically, as good as the un-rotated one: the baryon energy is degenerate in rotations. 

In quantum mechanics, however, the rotations are quantized. As first pointed out by Witten [2] and then derived 
using different techniques by a number of authors [10], the quantization rule is such that the lowest baryon multiplets 
are the octet with spin 1/2 and the decuplet with spin 3/2 (i.e. exactly those observed in nature) followed by the exotic 
aniz-decuplet with spin 1/2 again. The parity of all rotational states is the same. Those baryons are distinguished by 
the rotational wave functions depending on the eight Euler angles parameterizing the SU (3) rotation matrix R; the 
wave functions are given explicitly in Section IV. 

Qualitatively, one can think of different baryons as "living" in different parts of an 8-dimensional globe parametrized 
by 8 "Euler angles". The + lives near the North pole of that globe, at least in the academic limit of large N c . 
The average polar angle for the rotational state corresponding to the + vanishes as 1 / y/N c , see section IV. D. 
Therefore, in the limit N c — > oo one can approximate the rotation by small (kaon) fluctuations about the North pole. 
Mathematically, it comes to the Callan-Klebanov "bound-state approach to strangeness" [11] where one studies the 
linear response of a nucleon to a small-amplitude kaon perturbation, or the KN scattering, to see if there is a + 
resonance. In such approach the narrow Q + does not exist, at least in the Skyrme model for the KN scattering, 
unless one extends the parameters of the model [12]. The Skyrme model for the self-consistent chiral field is, however, 
not realistic, and it is unclear what lesson can one draw from the existence or non-existence of a resonance in this 
particular dynamical model. 

Even more important, it is exactly the situation where the large N c limit can hardly be trusted. In reality at N c = 3 
the 9 + rotational wave function is spread over the whole 8-dimensional globe and is far from the "North pole" . A 
quantum-mechanical model of the situation has been suggested by Cohen [13] and Pobylitsa [14]; the model can be 
solved numerically at any N c [15]. It turns out that the energy levels at N c = 3 differ radically from their positions at 
N c — ► oo. Given this experience, we shall treat the + rotational wave function exactly at N c = 3, see section IV. At 
the same time we shall neglect the fluctuations of the chiral field about its mean field value since these are suppressed 
additionally as are any generic loop corrections. 

In this approach, all low-energy properties of baryons from the ^8, , ^10, and ^10, multiplets (including 

e.g. parton distributions at low virtuality) follow from the shape of the mean chiral field in the common or 'classical' 
baryon; the difference and splitting between baryons from those multiplets arise exclusively from the difference in 
their rotational wave functions. This difference can be immediately translated into the quark wave functions of the 
individual baryons, both in the infinite momentum [16, 17] and the rest [18] frames. In Section III we present a 
compact general formalism how to find the 3-quark, 5-quark, 7-quark ... wave functions inside the octet, decuplet and 
antidecuplet baryons, which is further detailed in Sections V and VI. In Section VII we find the quark wave functions 
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of the 3Q components in the octet and decuplct baryons. In the non-relativistic limit (implying a weak mean field), 
we obtain the old SU(6) quark wave functions for the octet and decuplet baryons but with well-defined relativistic 
corrections. The 5Q wave functions in the ordinary and exotic baryons can be also found explicitly [17, 18], see 
Section VIII. 

In Sections IX-XI we develop a formalism how to compute observables related to the 3- and 5-quark Fock components 
of baryons, and apply it in Section XII to estimate the nucleon axial constant and the transition matrix element of 
the strange axial current between the 6 + and the nucleon: it gives an estimate of the 6 + — ► KN decay width. The 
latter turns out to be very small, 2-4 MeV, although with a large uncertainty discussed in Section XIII. 

The essence of QCD with its spontaneous breaking of chiral symmetry is that adding a low-energy pseudoscalar 
meson (or a QQ pair) to a baryon is equivalent to rotating the vacuum state along the Goldstonc valley, meaning 
no change of the physical state. In order to separate the true QQ pairs in a baryon from those in the vacuum, one 
has to consider baryons in the Infinite Momentum Frame (IMF). In this and only this frame the true QQ pairs in a 
baryon have an infinite momentum as contrasted to those in the vacuum, which have a finite momentum. Therefore, 
an accurate definition of what are the 3-, 5-,... Fock components of baryons can be made only in the IMF. It also 
has the advantage that the vector and axial currents with a finite momentum transfer do not create or annihilate 
quarks with infinite momenta. The baryon matrix elements are thus non-zero only between Fock components with 
equal number of quarks and antiquarks. 

Since the + has no 3Q component it means that one has to calculate the matrix clement between the 5Q component 
of the 6 + and the 5Q component of the nucleon. In principle, one has to add also the 7Q — ► 7Q, 9Q — ► 9Q... 
transitions, but we neglect them in the present paper. To control this approximation, we compute, using the same 
technique, the nucleon axial constant g^N). In the (very crude) non-relativistic 3Q approximation to nucleons, this 
constant is approximately 5/3 = 1.667; taking into account the 5Q component of the nucleon moves it to the value of 
1.36 being already not too far from the experimental value <?a(-V) — 1.27 [20]. It should be noted that the summation 
of the contributions of any number of QQ pairs in the Relativistic Mean Field Approximation to nucleons moves 
gA(N) quite close to the experimental value [21]. In the 5Q approximation to the + — ► KN transition, we obtain 
Su(@ — > KN) w 0.14 — 0.2 leading to the estimate Te ~ 2 — 4 MeV. In this estimate, we neglect the quark exchange 
contributions to the + — > KN transition, which are potentially capable of reducing further the width. Qualitatively, 
the axial constant of the + — > KN transition is small because it is analogous not to the large nucleon axial constant 
itself but to the change of this constant as one goes from the 3Q to the 5Q contribution. 



II. THE EFFECTIVE ACTION 



The effective action approximating QCD at low momenta describes "constituent" quarks with the momentum 
dependent dynamical mass M(p) interacting with the scalar (S) and pseudoscalar (II) fields such that X 2 + II 2 = 1 
at spatial infinity. The momentum dependence M{p) serves as a formfactor of the constituent quarks and provides 
the effective theory with the ultraviolet cutoff. Simultaneously, it makes the theory non-local. The action is [1] 

d 4 pd 4 P ' 



Seff= j-^t^(P) ^(2n) 4 S^ip-p')-^M{pj(^(p-p')+m(p-p') l5 ) v/MMJ V(p'), (2) 

where ip,^ are quark fields carrying color, flavor and Dirac bispinor indices. In the instanton model of the QCD 
vacuum from where this action has been originally derived the function M(p) is such that there is no real solution of 
the mass-shell equation p 2 = M(—p 2 ), therefore quarks are not observable as asymptotic states, - only their bound 
states. However, this is not the true confinement. Unfortunately, the instanton model's M(p) has a cut at p 2 = 
corresponding to massless gluons left in that model. In the true confining theory there should be no such cuts. 
Nevertheless, such M(p) creates some kind of a soft "bag" for quarks. Contrary to the naive bag picture which does 
not respect relativistic invariance, eq. (2) supports all general principles and sum rules for conserved quantities. 

The scalar, pseudoscalar [22], vector and axial [23] mesons follow from the correlation functions computed from 
eq. (2). The light-cone quark wave functions of the pion and of the photon have been found in Rcf. [24]; the 
electromagnetic pion radius has been computed in the original paper [1]. 

Turning to baryons, the mean S, II field (called chiral field for short in what follows) in the full non-local formulation 
(2) has been found by Broniowski, Colli and Ripka [25]. It sets an example how one has to proceed in the model 
calculations. However, to simplify the mathematics we shall use here a more standard approach: we shall replace the 
constituent quark mass by a constant M = M(0) and mimic the decreasing function M(p) by the UV Pauli-Villars 
cutoff [5]. 
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III. BARYON WAVE FUNCTION IN TERMS OF QUARK CREATION-ANNIHILATION OPERATORS 



Let a, oJ(p) and 6, 6^(p) be the annihilation-creation operators of quarks and antiquarks (respectively) of mass 
M, satisfying the usual anticommutator algebra {a(p)a^(p')} = {fo(p)&t(p')} = (27r) 3 S^ (p — p') and annihilating 
the vacuum state a,6|0>= 0, <0|a^,fet — 0. For quarks, the annihilation-creation operators carry, apart from the 
3-momentum p, also the color a, flavor / and spin a indices but we shall suppress them until they are explicitly 
needed. The Dirac sea is presented by the coherent exponent of the quark and antiquark creation operators [16], 



coherent exponent for QQ pairs = exp (^J (dp)(dp / ) a) (p) W(p, p') b^(p')^j |0>, 



(3) 



where (dp) — d 3 p/ (2ir) 3 and W(pi, P2) is the quark Green function at equal times in the background S, II fields [16, 
17] (see Fig. 2); we shall specify the function W below. In the mean field approximation the chiral field is replaced 
by the spherically-symmetric self-consistent field: 

7r(x) =n-rP(r), n = x/r, S(x) = E(r). (4) 

On the chiral circle (to which we restrict ourselves for simplicity) II — n ■ r smP(r), S(r) = cos P(r) where P(r) is 
the profile function of the self-consistent field. It is fairly approximated by [3, 8] 

P(r) =2atanf4^ °' > 



''() 



M ' 



(5) 



where M « 345 MeV is the dynamical quark mass at zero virtuality, known to fit numerous observables within the 
instanton mechanism of the spontaneous chiral symmetry breaking [1] . 

The self-consistent chiral field (5) creates a bound-state level for quarks, whose wave function ip\ eY satisfies the 
static Dirac equation with eigenenergy E\ ev [3, 26, 27]: 



V>lev(x) = 



-ie jk (<r-n) 4 fc j(r) 



ti + h M sin P- j{M cos P + E lcv ) = 0, 
j' + 2j/r - j M sin P - h(M cos P - E lev ) = 0, 



(6) 



where i = 1,2 is the spin and j = 1, 2 = u, d is the isospin index. In the non-relativistic limit (E\ cv « M) the L = 
upper component of the Dirac bispinor h(r) is large while the L=l lower component j(r) is small. Solving eq. (6) 
for the self-consistent field (5) one finds that 'valence' quarks are tightly bound (E\ cv = 200 MeV) but the lower 
component j(r) is still substantially smaller than the upper one h(r), see Figs. 3,4. 





FIG. 3: The space profile of the self-consistent chiral field FIG. 4: Bound-state quark upper s-wave component h(r) 



P(r) in light baryons. One unit on the horizontal axis is 
r = 0.8/M = 0.46 fm. 



(solid) and the lower p-wave component j(r) (dashed) in 
light baryons. The three valence quarks have the energy 
Ei cv = 200 MeV each. 



The valence quark part of the baryon wave function is given by the product of N c quark creation operators that fill 
in the discrete level [16]: 



valence quarks wave function = JJ / (dp) F(p) a f (p) 

l n J 



(7) 



lor=l ' 



F(p) = l(dp')J — [u(p) 70 iP lcv (p)(2n) 3 5(p-p')-W(p,p')v(p') l0 ipic V (-p')], (8) 
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where V'iev(p) is the Fourier transform of eq. (6). The second term in Eq. (8) is the contribution of the distorted Dirac 
sea to the one-quark wave function. u a (jp) and ^(p) are the plane- wave Dirac bispinors projecting to the positive 
and negative frequencies, respectively. In the standard basis they have the form 



Mp)= r^7„„ ' Mp)=\ V 2 r L T7 ' uu = l = -vv, (9) 





where e = +y p 2 + M 2 and s CT are two 2-component spinors normalized to unity, for example, 

81= (j), S2= (i)' ,T = 1 ' 2 - ( 10 ) 

The full baryon wave function is given by the product of the valence part (7) and the coherent exponent (3) 
describing the distorted Dirac sea. Symbolically one writes the baryon wave function in terms of the quark and 
antiquark creation operators [16]: 

B[a\bt]= J] /(dp)F(p) t(p) cxp('|(d P )(dp')at(p)^(p,p')6 t (p')) |0> . (11) 

color— 1 

At this point one has to recall that the saddle point at the self-consistent chiral field is degenerate in global 
translations and global SU(3) flavor rotations (1) (the SU(3) breaking by the strange mass can be treated as a 
perturbation later). Integrating over translations leads to the momentum conservation: the sum of all quarks and 
antiquarks momenta have to be equal to the baryon momentum. Integration over rotations R leads to the projection 
of the flavor state of all quarks and antiquarks onto the spin-flavor state B(R) describing a particular baryon from 

the (8, i + ) , (lO, f + ) or (TO, i + ) multiplet. 

Restoring color (a = 1,2,3), flavor (/ = 1,2,3), isospin (j = 1,2) and spin (er = 1,2) indices, the quark wave 
function inside a particular baryon B with spin projection k is given, in full glory, by [16, 17] 

tt fc (B) = [dRB* k (R)e a ^f[ /(d Pn )i?J;F-'" CT "(p n )a^ / ^(p n ) 

J 71=1 J 

• exp(y'(dp)(dpOa^ CT (p)4^(p,p')4/'6 t ^ V '(p'))|0> . (12) 

Acting on the vacuum state |0> the operators a> create three 'valence' quarks at the bound-state discrete level with 
the wave function F, while the a) operators in the exponent create any number of additional quark-antiquark 
pairs whose wave function is W. Eq. (12) is thus a full relativistic field-theoretic description of baryons, involving an 
infinite number of degrees of freedom. 

Note that the three 'valence' quarks are antisymmetric in color whereas the additional QQ pairs appear in color 
singlets. The spin-flavor quark structure of a particular baryon arises from projecting a general QQQ + nQQ state 
onto the quantum numbers of the baryon in question; this is achieved by means of integrating over all spin-flavor 
rotations R with the rotational wave function B^(R) unique for a given baryon. 

The third row of the matrix R? , f = 3, introduces strange quarks both at the valence level and in the sea; hence 
hyperons with explicit strangeness will, generally have valence s quarks, and non-strange baryons will contain ss pairs, 
even though only the u, d quarks are affected by the chiral field (4), which is reflected by the fact that the valence- level 
wave function F and the pair wave function W have not full SU(3) but only isospin indices j = 1, 2 = u, d. 

Eq. (12) encodes an enormous amount of information as it is the generating functional for the quark wave functions 
in all Fock components of baryons from the lowest multiplcts. Expanding the coherent exponent to the th , 1 st , 
2 nd ... order one reads off the 3-, 5-, 7-... quark wave functions of a particular baryon from the octet, decuplet or 
antidecuplet. All this information can be put in a compact form because the Relativistic Mean Field Approximation 
is being used. 

To make this powerful formula fully workable, we need to give explicit expressions for the baryon rotational wave 
functions B(R), the valence wave function F J ' CT (p) and the QQ wave function in a baryon Wjf a , (p, p'). 



IV. BARYON ROTATIONAL WAVE FUNCTIONS 



In general, baryon rotational states B{R) are given by the SU(3) Wigner finite-rotation matrices [28], and any 
particular projection can be obtained by a routine SU (3) Clebsch-Gordan technique. However, in order to see the 
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symmetries of the quark wave functions it is helpful to use explicit expressions for B(R), and integrate over the Haar 
measure in eq. (12) explicitly. 

We list below the rotational D-functions for the multiplets (8, i), (10, |) and (10, |) in terms of the product of 
the R matrices. Since the projecting onto a specific baryon in eq. (12) involves its conjugate rotational wave function, 
we list the conjugate functions only. The un-conjugate ones are obtained by hermitian conjugation. 



A. (8,1) 

From the SU (3) group point of view, the octet of baryons transforms exactly as an octet of mesons; therefore, its 
rotational wave function can be composed of a quark (transforming as R) and an antiquark (transforming as R^). 
Accordingly, the rotational wave function of an octet baryon labeled by a = 1 ... 8 and having a spin index k = 1, 2 is 



D^*(R)] a ~ e kl R\ l (t a ) f g R g 3 . 



(13) 



where e k i is the antisymmetric 2x2 tensor and t a are the SU(3) generators. In particular, the proton (a = 6 + i7) 
and neutron (a = 4 + i5) rotational wave functions with spin k = 1, 2 are 



p k (R)* = V8e kl R\ l Rl n k (R)* = V8e kl rI 1 R%. 



(14) 



B. (10,|) 

The decuplet states can be composed of three quarks; they are labeled by a triple flavor index {/1/2/3} symmetrized 
in flavor and by a triple spin index {fcife^s} symmetrized in spin: 



£(10,§)* (i?) ' 



{/l/2/ 3 },{felfc2fc 3 } 

For example, the A-resonance rotational wave functions are 

3 



e k ' lkl e k < 2k 2 e k< 3k3 -ft^ -K/ 2 Kf 3 



sym in {/1/2/3} 



A ++ , spin projection + 
A , spin projection + 



a++(r)* = VTor\ 2 r\ 2 r\ 2 , 



A?(i?)' 



10^ 2 (2i?l 2 i?^ 1 +i?| 2 i? t 1 1 ). 



(15) 

(16) 
(17) 



10, 



From the SU(3) group point of view, the antidecuplet can be composed of three antiquarks and its conjugate 
rotational wave function is 



D^*(R) 



{/1/2/3} 



p/i 0/2 73/3 
K 3 H 3 H k 



sym in {/1/2/3} 



In particular, 



Q + , spin projection k 
neutron* from 10 , spin projection k 



Q k (R)* = VZOR^RlR 3 
nf(RY : 



k > 
2 d3 d3 



10 (2R 3 R 3 R k + R 3 R 3 R k ) . 



(18) 



(19) 
(20) 



All the rotational wave functions above are normalized in such a way that for any (but the same) spin projection 

r dRB: pin (R)B s ^(R) = l: (21) 



for different spin projections the integral is zero. Rotational wave functions belonging to different baryons are also 
orthogonal. It can be easily checked directly using the concrete parametrization of the SU(3) rotation matrices R 
from Appendix A and performing the 8-dimensional integration with the measure defined there. 



7 



D. Large N c limit 

If N c is not equal to three but is treated as a free parameter, the lightest baryons are not the octet, decuplet 
and antidecuplet but some large ST/ (3) multiplets whose dimensions depend on N c . What SC/(3) multiplets are the 
large- A c prototypes of the usual multiplets at A c = 3, is not uniquely defined. It seems natural to define the prototype 
multiplets in such a way that their lightest members are "nucleons" with spin and isospin |, "A's" with spin and 
isospin |, and "O" with spin | and isospin 0: this prescription is sufficient to define unambiguously the large- N c 
prototypes of the octet, decuplet and antidecuplet [13, 29, 30]. 

The rotational wave functions of the large- A c analogs of the N, A and 6 are obtained from eqs.(14,16) and (19) 

by multiplying the corresponding equations by a factor (-R 3 ) " 3 . In Appendix A we give a concrete example of the 
parametrization of a general SU(3) rotation matrix R in terms of eight "Euler" angles. In fact they parameterize the 
S 3 x S 5 space, - the direct product of the 3d and 5d spheres. In this parametrization, 

Rl = e ia21 cos & cos 0, 0, <j) 2 S (0, |) , (22) 

where and (f>2 can be viewed as polar angles of the 5d sphere. It is clear that at N c — > 00 the rotational wave functions 
of the "A" , "A" and "9" are squeezed near the "North pole" of the sphere S 5 since the average polar angles vanish 
as 0,(f>2 ~ 1/s/Nl. The rotated self-consistent field (1) can be also parameterized a la Callan-Klebanov [11]: 

U = RVrf = Wu K W, V(x) = ( CXP [z(n ' T)P(r)] ) , (23) 

1 

where the meson SU(3) unitary matrix Uk is, for small meson fluctuations <fi about the self-consistent field V, 

U K = U + i(t> A X A , A = \..&, (24) 

One can compare both sides of eq. (23) and find the meson fields in baryons corresponding to rotations. In particular, 
for rotations "near the North pole" i.e. at small angles 0, <p2, one finds the kaon field 

P(r) P(r) 
K+ = -V2sm^-[0n z + (f>2{n x -iny)}, K° = —V2 sin —j^- [0(n x + in y ) — (f>2n z ] , (26) 

meaning that at large N c the amplitude of the kaon fluctuations in the prototype baryons "A", "A" and "0" is 
vanishing as ~ Therefore, the Q problem becomes that of the linear response of a nucleon to a small kaon 

fluctuation, and can be studied in a particular model for the effective chiral lagrangian [12]. However, in reality at 
A c = 3 the rotational wave functions of A (14), A (16) and 6 (19) correspond to large angles 0, <p 2 and arc not 
concentrated near the "North pole" . It means that the kaon field in these baryons is generally not small. Therefore, 
in what follows we shall use the exact /V c = 3 rotational wave functions (14,16,19). 



V. QQ PAIR WAVE FUNCTION 

As explained in Refs. [16, 17], the pair wave function Wj? a , (p, p') is expressed through the finite-time quark Green 
function at equal times in the external static chiral field (4); schematically, it is shown in Fig. 2. We shall need the 
Fourier transforms of the self-consistent chiral field, 

II(q)j, = J d 3 xe- iq x (n • t)], sin P(r), £(q$ = J d 3 x e - lq ' x (cosP(r) - l)<5j,, (27) 

where IT(q) is purely imaginary and odd while S(q) is real and even. 

In Refs. [16, 17] a simplified interpolating approximation for the pair wave function W has been derived, which 
becomes exact in three limiting cases: i) small pion field P(r), ii) slowly varying P(r), iii) fast varying P(r). In the 
infinite momentum frame the result is a function of the fractions of the baryon's longitudinal momenta carried by the 
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quark (z) and antiquark (V) of the pair, and their transverse momenta p_L,p'j_ [31]: 

MM 



7A^P±;z\p' ± ) = {xi,(q)[M(z'-z)T S + (zp'-z'p) ± .T ± Z, 

+ i n^,(q)[-M(z' + z)l + ie afi {zp' -z'p) ia ri£,}, (28) 
Z = M 2 zz'(z + z') + z(p'l + M 2 ) + z'ipl + M 2 ), q = ((p + p')x, + z')M) . 
Here j,j' — 1,2 are the isospin and a, a' = 1,2 are the spin projections, t\^,3 & re Pauli matrices, e aj a is the 



mq),Iiq) 
q=p+p' 




p,j, a 
quark 



antiquark 
P'J'. * 



FIG. 5: The QQ pair wave function in a baryon in the Relativistic Mean Field Approximation is related to the Fourier transform 
of the static self-consistent chiral field. 

antisymmetric 2d tensor, the primed indices refer to the antiquark; M is the baryon and M is the constituent quark 
masses, q is the 3-momentum of the pair as a whole, transferred from the background field S(q), n(q). 

The pair wave function W is normalized in such a way that the creation-annihilation operators in eq. (12) satisfy 
the anticommutation relations 

K^(z l!m ),4 2/2CT > 2 ,p 2 x)} - S%S%S% S( Zl z 2 ) {2n) 2 5^ 2 \p 11 _ p 2X ) (29) 

and similarly for b, b', and the integrals over momenta there are understood as fdzj d 2 p±/ (2ir) 2 . 

The pair wave function can be written in a more compact form by introducing the fraction of the longitudinal 
momentum of the pair carried by the antiquark y 7 and the transverse combination Q^, 

z' z + z' zp\ - z'p_L 
V=^— n Qz = — r-j-, Q± = —j • 30) 

z + z' M z + z' 

With this substitution eq. (28) takes the form 

W^AV^M = — n Qi + M* + y(l-y)<? " (31) 



VI. BOUND-STATE WAVE FUNCTION 

As seen from eq. (8), the discrete-level wave function F J<T (p) = F^(p) + F^ a (p) consists of two pieces: one is 
directly the wave function of the valence level, the other is related to the change of the number of quarks at the 
discrete level as due to the presence of the Dirac sea; it is a relativistic effect and can be ignored in the non-relativistic 
limit, together with the lower L = 1 component j(r) of the level wave function. Indeed, in the baryon rest frame the 
evaluation of the first term in eq. (8) gives 



™=^v^W^<4 (32) 

where h(p),j(p) are the Fourier transforms of the valence wave functions (6): 

h{p) = J d 3 xe- ip -*h(r) =4tt J 'drr 2 ^^h{r), (33) 

f p a An [ 

j a (p) = J d 3 xe~ lp x (-in a )j(r) = J^jip), j(p) = ~^ J dr (P r cos P r - sin pr)j(r). (34) 
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One sees that the second term in eq. (32) is double-suppressed in the non-relativistic limit E\ cv w M: first, owing to 
the kinematical factor, second, since in this limit the L — l wave j{r) is much less than the L = wave h(r). 
In the infinite momentum frame the evaluation of the bispinors u, v from eq. (9) produces [16, 17] 



e^h(p) + (p z l + ie aP p±aT±p) a al e jV 3 M 



p z =zM-Ei ev 

Similarly, the evaluation of the "sea" part of the discrete-level wave function gives 



(35) 



(36) 



p z =zM-E lo 



where the pair wave function (28) has to be used. 

In the following evaluation of the baryon matrix elements we shall neglect the relativistic effects in the discrete level 
wave function replacing it by the first term in eq. (35): 



e j °h(p)\ _ 



p z =zM-E l , 



(37) 



We have now fully determined all quantities entering the master eq. (12) for the 3,5,7... Fock components of baryons' 
wave functions. 



VII. 3-QUARK COMPONENTS OF BARYONS 

The absolute majority of baryon models focus on the 3-quark Fock components of the usual (non-exotic) baryons. 
We have already mentioned in the Introduction that it is a crude approximation to reality: the 5-, 7-,... quark 
components in the nucleon are not only non-negligible but critical for explaining such important characteristics as the 
nuclcon a term or the fraction of nucleon spin carried by quarks. Nevertheless, the 3-quark component is definitely 
important. In this section we derive the 3Q wave functions of the octet and decuplet baryons from our master 
equation (12) and show that in the non-relativistic limit they become the well-known SU(6) wave functions of the old 
constituent quark model. 

One gets the 3Q component of a baryon by ignoring the coherent exponent with QQ pairs in eq. (12); each of the 
three valence quarks is rotated by the matrix Rj where / = 1, 2, 3 is the flavor and j = 1, 2 is the isospin index. To 

obtain the color-flavor-spin-space 3Q wave function of a particular baryon from the ^8, ^ + ^jor the ^10, | one has 

to integrate in eq. (12) over all 8-parameter SU(3) rotations R with the (conjugate) rotational wave function B£(R) 
corresponding to the chosen baryon with spin projection k. These functions are given in Section IV. The arising 
SU(3) group integrals are of the type 

T (B)^t k = jdRBt{R)R^R%R% (38) 

where the three unitary matrices , R^ , R^ rotate the flavor of the quarks on the discrete level. These tensors are 
computed in Appendix B: for baryons from the ^8, -| + ^thc relevant integral is eq. (B8), and for the ^10, | + ^jit is 
eq. (B12). The tensor T must be now contracted with the three discrete-level wave functions from Section VI 

( Pl )F^ (p 2 )F^ (p 3 ). (39) 

In general the 3Q wave function depends on all four quark "coordinates": the position in space (r) (or the three- 
momentum p), the color (a), the flavor (/) and the spin (<r), and also on the baryon spin projection k. The wave 
function must be antisymmetric under permutation of all four "coordinates" for a pair of quarks. We suppress the 
trivial color wave function e aill2a3 which factors out. In the non-relativistic approximation we use the simplified level 
wave function (37) and for clarity pass back to the coordinate space. We thus obtain, for example, the 3Q wave 
function of the neutron: 

(|n> fe ) /l/2/3 ' CTl<T2CT3 ( ri ,r 2 ,r 3 ) - ||e /l/2 S f 2 s 6? h(n)h(r 2 )h(r 3 ) 

+ permutations of 1,2,3, (40) 
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times the antisymmetric e QlQ2a3 in color. In this equation the flavor indices assume only two values: /i,2,3 = 1, 2 = u, d. 
Eq. (40) says that the neutron spin is carried by the d-quark, and the ud pair is in the spin- and isospin-zero 
combination. It is better known in the form 

K> = 2d| (ri)dT (r 2 )t4 (r 3 )-dT (ri)«T (r 2 )<U (r 3 )-«T (ri)<U (r 2 )dT (r 3 ) 

+ permutations of n, r 2 , r 3 , (41) 

which is the well-known non-relativistic SU(6) wave function of the nucleon, with a concrete space distribution of 
quarks, shown in Fig. 4. 

Similarly, the 3Q wave function of the A resonance with spin projection 1/2, which may be compared with 
that of the neutron, can obtained from the group integral (B12), and reads 



|A°T> /l/2/3 ' <7l,T2CT3 (ri,r 2 ,r 3 ) = ^ (sfrM* + S^SfS? + SM'ti") 

■ (5? 6?%* + 8? 8? 8? + ^C^ 3 ) /i(n)/i(r 2 )/i(r 3 ) (42) 

which can be also presented as a familiar SU (6) wave function 

| A 1> = «T (ri) d| (r 2 ) dj. (r 3 ) + dj. (n) u| (r 2 ) d? (r 3 ) + d] (n) d| (r 2 ) uj. (r 3 ) 

+ permutations of n, r 2 , r 3 . (43) 

There are, of course, relativistic corrections to these S'C/(6)-symmetric formulae, arising from i) exact treatment of 
the discrete level, eqs.(35,36), and ii) additional QQ pairs described by eq. (28). Both effects are generally not small. 

VIII. 5-QUARK COMPONENTS OF BARYONS 

One gets the wave functions of the 5Q component of baryons by expanding the coherent exponent in the generating 
functional (12) to the linear order in the QQ pair. The SU(3) group integral involves now three i?'s from the level 
and R, W from the pair, times the (conjugate) rotational wave function B^{R) of the baryon in question: 

T (B)^±H k = J dR B k ( R ) R h R h R h R u R \T ■ ( 44 ) 

We shall systematically attribute the indices 1,2,3 to the valence quarks, index 4 to the extra quark of the QQ pair, 
and index 5 to the antiquark. The group integral (44) is computed in Appendix B: for octet baryons the result is 
given in eq. (B14) and for the antidecuplet baryons it is given in eq. (B17). To obtain the 5Q wave function of a 
baryon, one has to contract T from eq. (44) with three valence quark wave functions F (39) and with the pair wave 
function W (28). 

In general, the 5Q wave functions look rather complicated as they depend on five quark "coordinates" , including 
their coordinates proper (or 3-momenta), spin, flavor and color. We do not write explicitly the color degrees of freedom 
but always imply that the (1, 2, 3) quarks of the level are antisymmetric in color while the quark-antiquark pair (4, 5) 
is a color singlet, as it follows from eq. (12). For example, the 5Q wave function of the neutron is 

(|n> fc )^ 5 /3/ ^ 2CT ^(pi---P5) - ^^ 1CT1 (pi)^ 2a2 (p2)^ 3CT3 (P3)v^;;;(p4,P5) 

+ e f ^e jlj4 [4 2 ^ (4$6&-8£6£) 

+ permutations of (1,2,3). (45) 



Terms of the type of 5j? mean the flavor-symmetric combination ss + uu + dd, however quarks from this combination 
are partly inside the pair wave function W but partly in the "valence" bound state. We have not invented how to 
present it in a more compact form; however, eq. (45) is a working formula which allows to get compact physical 
results, see Section XII. The 5Q wave function of the proton is the same, with the replacement 8 2 1 ' 2 ' iA — > <5^ 1,2 ' 3,4 , 
meaning that one d-quark must be replaced by the u-quark. 
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Turning to the exotic baryons from the (^10,^ J, projecting the three quarks from the discreet level onto the 

antidecuplet rotational function (18) gives an identical zero in accordance with the fact that the exotic antiducuplet 
cannot be made of 3 quarks, see eq. (B16). The non-zero projection is achieved when one expands the coherent 
exponent at least to the linear order. For example, one gets then from eqs.(19,B19) the 5Q wave function of the 
9+: 



\&+>^ f ^^ (Pi...p 5 ) = ^W3X^ 2e , 3 ^^ 

+ permutations of (1, 2, 3). (46) 

The color structure of the antidecuplet wave function is e ai 012013 5^. The quark flavor indices are fi^ = 1,2 = u, d, 
and the antiquark is s owing to 5^ . Naturally, we have obtained the quark content 9 + = uudds where the two (ud) 
pairs are in the isospin-zero combination, thanks to e-' 1 * 2 e-' 3 -' 4 . 

To make contact with other work where the + wave functions were obtained in various non-relativistic models 
or discussed in that framework [32], one has to pass to the coordinate space and write eq. (46) in the 6 + rest frame 
using the non-relativistic approximation (37) for the level wave function. We obtain 



>/5^ /3/4 ^ ff2CT3ff4 (n--- r 5) = ^e^ef^Sle^h( ri )h(r 2 )h(r 3 )W^(v 4 ,r 5 ) 

+ permutations of (1,2,3) (47) 

where the pair wave function in the coordinate space W / (r 4 ,r 5 ) can be found in Ref. [18]. The structure e^ 2 e CTl<T2 
clearly shows that there is a pair of ud quarks in the spin and isospin zero combination, exactly as in the nucleon, 
eq. (40). However, it does not mean that there are prominent scalar isoscalar diquarks either in the nucleon or in the 
+ : that would require their spatial correlation which, as we see, is absent in the mean field approximation. The QQ 
pair wave function W is a combination of four partial waves with different permutation symmetries, in accordance 
with the general analysis by Bijker, Giannini and Santopinto, Ref. [32]. The amplitudes of those partial waves depend 
separately on the coordinates r 4 5 measured from the baryon center of mass. More explicit formulae are given in 
Ref. [18].' 



IX. THREE QUARKS: NORMALIZATION, VECTOR AND AXIAL CHARGES 

The normalization of a baryon wave function in the second-quantization representation (12) is found from 

Af(B) = lj, fc <*t BZ *f > . (48) 

The annihilation operators in must be dragged to the right where they ultimately nullify the vacuum state |0> 

and the creation operators from should be dragged to the left where they ultimately nullify the vacuum state <0|. 
The result is non-zero owing to the anticommutation relations (29) or the "contractions" of the operators. 

For the 3Q Fock component of a baryon, there are 3! possible (and equivalent) contractions, and the ensuing 
contraction in color indices gives another factor of 3! = e aia2a3 e aia2 a 3 . Flavor projecting to a baryon with specific 
quantum numbers gives the tensor (38), or its hermitian conjugate for the conjugate wave function. Hence the 
normalization of the 3Q component, shown schematically in Fig. 6, left, is 

■(2n) 2 5((p 1+ p 2 +p 3 ) ± ) F^{ Pl )F^(p 2 )F^{p 3 ) Fl^Fl^F^ips) (49) 

where F JCF (z,p±) are the level wave functions (35,36). In the non-relativistic limit F ja (p)Fi cr (p) ~ Sf h 2 (p), see 
eq. (37). Therefore in this simple case the normalization is the full contraction of the two T tensors, times an integral 
over momenta which can be performed numerically once the level wave function h(p) is known. 

A typical physical observable is the matrix element of some operator (which should be written down in terms of 
the quark annihilation-creation operators a, b, a\tf) sandwiched between the initial and final baryon wave functions 
(12). We shall consider as examples the operators of the vector and axial charges which can be written through the 
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Q 



annihilation-creation operators as 

b^ h p(z, P± )b e7I (z,p^,r h 



dz 



(2^) 2 



(-0-3); 



("OS)? 



(50) 



where is the flavor content of the charge, and it, p = 1,2 are helicity states. For example, if we consider the p + = du 
current which annihilates u quarks and creates d quarks and annihilates d antiquarks and creates u quarks, the flavor 
currents in eq. (50) are Jf l (p + ) — <5|<5^. Notice that there are no a)tf or ab terms in the charges. This is a great 
advantage of the infinite momentum frame where the number of QQ pairs is not changed by the current. Hence there 
will be only diagonal transitions between Fock components with equal numbers of pairs, see Fig. 6, right. 




FIG. 6: Graphs showing the normalization of a 3-quark 
component of a baryon (left) and the matrix element of 
a local operator denoted by a circle (right). 



FIG. 7: Direct (left) and exchange (right) contri- 
butions to the normalization of the 5-quark com- 
ponent of a baryon. The upper rectangles denote 
QQ pairs. 



In the matrix elements between the 3Q components the b^b part of the current is passive as there are no antiquarks. 
The a^a part is a sum over colors. As in the normalization, one gets the factor 6 • 6 from all contractions. Let it be 
the third quark whose charge is measured: there is a factor of 3 from three quarks to which the charge operator can 
be applied, see Fig. 6. Denoting for short J(dpis) the integrals over momenta with the conservation <5-functions as 
in eq. (49) we arrive at the following expression for the matrix element of the vector charge: 



2) = 



(6-6-3) 



[F jiai ( Pl )F j 



(dpi-a) 



2) F i 3T3 (pa) 



XT3 T93 



(51) 



One can easily check using eq. (14) that, say, for the p — > np + transition, the above vector charge gives exactly the 
same expression as for the normalization (49). Therefore, the gy of this transition is unity, as it should be for the 
conserved vector current. 

We consider here for simplicity only matrix elements of operators with zero momentum transfer. If it is non-zero, 
the generalization is obvious: one has to change the momentum of one of the quarks on which the operator acts, by 
the corresponding momentum transfer, and leave the rest quarks momenta unaltered. 

For the axial transition, one replaces averaging over baryon spin by |(— 03); , and the axial charge operator is now 
(— (Ta)^ 3 instead of 5^, see eq. (50). All the rest is the same as in eq. (51): 



2) = 



(6-6-3) 



^(--3)?T(l)£j^ fc T(2)^ 



(dp 



1-3) 



[ F h<n (p^F^ (p 2 )F^ (pa)] F} iai {px)F} 2(T2 ( P2 )Fl T3 (P3) J 93 



(52) 



The result, however, is now different as the axial charge is not conserved. For example, for the p — > mr + transition one 
gets the expression identical to that for the normalization but with the factor 5/3. It means that we have obtained 

in the non-relativistic limit for the 3Q component of the nucleon g^(N) — 5/3. It is the well-known result of the 
non-relativistic quark model. However, it is modified by the relativistic corrections to the valence quark wave functions 
(35,36) and by the 5Q component of the nucleon. 
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X. FIVE QUARKS: NORMALIZATION, VECTOR AND AXIAL CHARGES 



Already in the normalization of the 5Q Fock component of a baryon there are two types of contributions: direct 
and exchange ones, see Fig. 7. In the former, one contracts from the pair wave function with an a in the conjugate 
pair, and all the valence operators are contracted with each other. There are 6 such possibilities, and the contraction 
in color gives a factor 3 • 6, all in all 108. In the exchange contributions, one contracts cv from the pair with one 
of the three a's from the valence level. Further on, a from the conjugate pair is contracted with one of the three 
a^'s from the valence level. There are 18 such possibilities but the contraction in color gives now only a factor of 6. 
Therefore for the exchange contractions we also get a factor of 108 but with an overall negative sign as one has to 
anticommute fermion operators to get the exchange terms. As a result we obtain the following general expression for 
the normalization of the 5Q Fock component: 



F h<n ( Pl )F^ ( P2 )F^ (p 3 ) W]f a l ( P4 , P5 ) F? iai ( Pl )Fl a2 (pa) 
>,U (Ps) W&(P4, Pb) S% 5% - F} 3a M) W ofJ 3 (P3 , P^fA 



where we have denoted 

'(dpi-:, 



J(dpi- 5 ) = Jdz 1 _ 5 5(1 - zi - ... - z 5 ) J d ^ro^ ( 2?r ) 2 S (P^ + - + P6-i-)- 
The flavor tensor here is the group integral projecting the 5Q state onto a particular baryon, see eq. (44). 



(53) 



(54) 



X 



I IXi 



FIG. 8: Direct contributions to the matrix ele- FIG. 9: Four types of exchange contributions to the matrix 
ment of an operator, in the 5-quark component of element in the 5-quark component of a baryon. 
a baryon. The operator is applied to the antiquark 
(left), to the quark from the pair (middle) and to 
the quark from the valence level (right). 

The ratio of the normalization factors AT^ /Af^ gives the probability to find a 5Q component in a mainly 3Q 
baryon. It depends on the mean field inside a baryon through the pair wave function W (and is quadratic in the mean 
field), and on the particular baryon through its spin- flavor content T. 

For the vector and axial transitions there are three basic contributions: one when the charge of the antiquark is 
measured, the second when the charge operator acts on the quark from the pair, and the third when it acts on one of 
the three valence quarks. These three types are further divided into the direct and exchange contributions (Figs. 8,9). 
We write below only the direct contributions; the exchange ones can be easily constructed from the graphs in Fig. 9. 

The vector transition: 



y(5)direct/j _^ 2) 



108 



2 6 i 1 ^hxhhuJ^k 1 Whhgag*M J 1*1-5) 
pirn ( pl ) F h* 3 ( P2 ) F ^ (pg) (p 4 , pb) F/ icti (pi)Fl^ ( P2 )Fl T3 (p 3 ) (p 4 , ps) 



(55) 



The axial transition: 

^(5)direet/j _^ 2) 



108 



(dpi -5 



( - \k rpf-i \flhf3fi,]5 rpfr)\lll2hU,g5,l 

2 I a 3)l 1 \ L )jij a j 3 jA,f B ,k 1 \ z )fxf2g39i,h 
F h*r fa)F*°* ( P2 )F^ (p 3 ) (p 4 ,ps) ( Pl )F/ 2CT2 (JP2)F} 3T3 (p 3 ) W l c \ T J A (p 4l p 5 ) 



(56) 
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where is the flavor content of the current defined in the previous section. 

In the next sections we apply these general formulae to the calculation of the nucleon axial charge and the 6 + 
width. 



XI. FIVE QUARKS: OVERLAP INTEGRALS IN THE INFINITE MOMENTUM FRAME 



It takes a few minutes by Mathematica to perform the contractions in eqs. (53,55,56) over all flavor (/, g), isospin 
(j, I) and spin (a, r) indices. After all contractions are performed, one is left with scalar integrals over longitudinal (z) 
and transverse (p_i_) momenta of the five quarks. The integrals over the relative transverse momenta in the QQ pair 
are generally UV divergent, reflecting the divergence of the negative-energy Dirac sea of quarks (Fig. 1). In reality, 
this divergence is cut by the momentum-dependent dynamical quark mass M(p), see eq. (2). Following Rcf. [5] where 
parton distributions in nucleons have been computed, satisfying all general sum rules, we shall mimic the fall-off of 
M(p) by the Pauli-Villars cutoff at Mpy = 557 MeV (this value is chosen from the requirement that the constant 
F-K = 93 MeV is reproduced from M(0) = 345 MeV). 

The pair wave function W (28) is determined by the Fourier transforms of the mean chiral field II(q) and £(q) 
(27): we find 



n(q)i 



n(g) 



4tt r 



dr sin P(r)(qr cos qr — sin qr) < 0, 



47T r°° 

E(q)5 = 5) £(«), £(?) = — / drr (cos P(r) - 1) sin qr < 0. 

\Q\ Jo 



(57) 
(58) 



Actually q is the 3-momentum of the QQ pair, which in the infinite momentum frame is q = (p4_L +Ps_l, (z& + zc,)M.). 

In the "direct" 5Q — ► 5Q transitions (53,55,56) with zero momentum transfer the following four scalar integrals 
arise from squaring eq. (31), corresponding to i) the full square of II(q), ii) the square of X(q), iii) the square of the 
third component IL3(q), and iv) the mixed IT3(q)£(q) term: 



K nn = 



M2 
2tt 

2tt. ; 



_d 3 q 

(2- 



" d 3 q ( q z 
(2^ $ U' Q 



) e(q z ) qz X 2 (q) f dy 



M_ 2 



d 3 q 



(2tt) 3 \M : 
d 3 



% z )^n^(q) dy 



d 2 Q± 
(2tt)2 

d 2 Q± 
(2n) 2 

d 2 Q A 

(2ny 



Ql+M 2 



(Ql + M 2 +y(l-yWY 
Ql + M 2 (2y-1) 2 

(Qi+M2 + y(l_y) q 2) 2 

Ql + M 2 



(M- 



- (M- 



>M PV ) 
>M PV ) 



d 2 Q A 



(Ql + M 2 + y{l-y)q 2 Y 
Ql + M 2 (2y-1) 



(27T) 2 



(Ql + M 2 +y(l-y)q 2 ) 



(M^Mpy) 
, -(M^Mpy) 



(59) 
(60) 
(61) 
.(62) 



We have rearranged the integrals dpi-5 such that we first integrate over the relative momenta inside the QQ pair 
y, Q± (see eq. (30)) and then over the 3-momentum q of the pair as a whole. As explained above, we regularize all 
integrals over the relative momenta by the Pauli-Villars subtraction. The step function 0(q z ) ensures that in the IMF 
the longitudinal momentum carried by the pair is positive. By <j>(z,qj_) we denote the probability that three valence 
quarks "leave" the longitudinal fraction z = Z4 + z 5 = q z /M and the transverse momentum q^ = + p 5 j^ to the 
QQ pair: 



$(z,qj-) 



d 2 P 



1,2,3 -L 



dzi,2,3 (2tt) S(p 1± + p 2 ± + p 3 ± + q±) 6(1 - 
Vz 1,2,3 = zx.2, 3 M - E lcv . 



(2tt) 6 

Pi, 2, 3 = (p± 1,2,3, Pz 1,2,3), 



zx-Z2- z 3 ) h 2 ( Pl )h 2 (p 2 )h 2 (p 3 ) 1 (63) 



In the 3Q components of baryons, there are no additional QQ pairs, and all quantities considered in Section IX are 
proportional to $(0,0). Since the normalization of the discrete-level wave function h(p) is arbitrary, we choose it 
such that $(0,0) = 1. 



Let us give a few examples how the normalization, vector and axial charges of the 5Q components of baryons are 
expressed through the integrals (59-62) after all contractions in eqs. (53,55,56) are performed. 
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Nucleon normalization: 

1 8 

tf 6) (N) = — (lltf™ + 2AK aa ). (64) 
5 

For the vector charge of the n — > p transition one gets exactly the same expression, which demonstrates that the 
vector charge is conserved in each Fock component separately. The vector charge of the 9 + — > K+n transition turns 
out to be identically zero: it reflects the known fact that matrix elements of SU(3) flavor generators between different 
irreducible representations, in this case between 10 and 8, are zero; it serves as an additional check of eq. (55) since 
individual contributions in that equation are non-zero. 



Nucleon axial charge: 

A^(p -» 7r+n) = 4(209if™ + 559K aa - 34K 33 - 356K 3a ). (65) 
25 

+ normalization: 

N^(Q) = —{K me + K aa ). (66) 
5 

Axial charge of the 6 + — ► K+n transition: 

A^(B+ - K+n) = ^^(-7^ - hK aa + 8K 33 + 28X 3ct ). (67) 

Notice that the coefficients are an order of magnitude less in the + than in the nucleon case, ft should be noted 
that we have independently derived eqs. (64-67) in another way by applying the charge operators directly to the five 
quarks and using the SU(3) Clebsch-Gordan machinery for projecting the 5Q states onto the baryons in question. 
Since this technique is different from the one presented here, it serves as a powerful check of the above expressions. 
We now proceed to evaluate them numerically. 



XII. NUMERICAL RESULTS 



For the numerical evaluation of the integrals involved in the 5Q matrix elements we use the quark mass M = 
345 MeV, the self-consistent profile function (5), the Pauli-Villars mass Mpy = 557 MeV, and the baryon mass A4 = 
1207MeV, as it follows for the "classical" mass (i.e. without quantum corrections) in the mean field approximation 
[8]. The self-consistent pseudoscalar II(q) and scalar £(q) fields, as given by eqs. (57,58) are plotted in Fig. 10. The 
probability distribution &(z,c[±) (63) that the QQ pair carries the fraction z of the baryon momentum and the 
transverse momentum is plotted in Fig. 11. 

With these functions, the numerical evaluation of the integrals (59-62) yields 

K nv = 0.0623, K aa = 0.0284, K 33 = 0.03 72, K 3a = 0.0333. (68) 
Putting these values into eqs. (64-67) we obtain 
Nucleon 5Q normalization: 

Af^(N) = 4.813. (69) 

Nucleon 5Q axial charge: 

A (5) (p^7r+n) = 3.779. (70) 

Q + 5Q normalization: 

7V (5) (9) = 0.652. (71) 

0+ 5Q axial charge: 

A (5) (@+ ^ K + n ^ = o 607 (72) 
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FIG. 10: The self-consistent pseudoscalar — |q|n(q) 
(solid) and scalar — |q|E(q) (dashed) fields in baryons in 
the momentum representation. The horizontal axis is in 
units of M = 345 MeV. 




FIG. 11: The probability distribution <E>(z,qx) that the 
QQ pair carries the fraction z of the baryon momentum 
and the transverse momentum q^ plotted in units of M = 
345 MeV. 



One has to add the 3Q nucleon normalization computed from eq. (49) 

A/" (3) (7V) = 9$(0,0) = 9 (73) 
and the 3Q nucleon axial charge computed from eq. (52) 

A^(p^ n+n) = 15$(0,0) = 15, (74) 
from where it follows that in the non-relativistic 3Q approximation the nucleon axial charge is 

(3) A^jp^ir+n) 5 

9a{n) = jmw) = 3 * L67 (75) 

which is the well-known result of the non-relativistic quark model. 
In the 5Q approximation, the nucleon axial charge is 

(S), m ^ (p -> n+n) + (p -> Tr+n) oa 

9a {n) = — jmiN^mw) — m 

which brings it closer to the experimental value gA(N) = 1.27. The account for any number of pairs and for relativistic 
corrections in the 1/N C expansion brings qa very close to the experimental value [21]. 
We note that the ratio of the 5Q to the 3Q normalization in the nucleon is 

0.535 .50%. (77) 

On the one hand, it means that the 5Q Fock component of the nucleon is quite substantial but on the other hand it 
implies that antiquarks carry roughly only 

-^-j 7% 

of the nucleon momentum, assuming the antiquark carries 1/5 of the momentum in the 5Q component [33]. We have 
not evaluated the 7Q... normalization in the nucleon (which would follow from expanding the coherent exponent in 
eq. (12) to higher orders) but expect that higher Fock components are suppressed, roughly, by factorials following 
from the expansion of the exponent. At large iV c , however, there would be on the average 0(N C ) QQ pairs in the 
nucleon. 



17 



Turning to the axial constant of the 6 + — > KN transition we obtain 



4(5) (pi+ . ir+ rl \ 

g A {Q-*KN)= - \ ' = 0.202 (78) 

y/NW(e) V^ {3) (N) +7Ffi{N) 

being substantially less than the nucleon axial charge computed in the same approximation. The quantity is similar 
in spirit (and magnitude) not to the nucleon axial coupling itself but to its change as due to the 5Q component in the 

nucleon, g^\N) — g^p(N) — 0.31. It is additionally suppressed by the S£/(3) group factors for the 10 — > 8 transition. 

Assuming the approximate ST/ (3) chiral symmetry (which was the base for using the + wave function (46) in the 
first place) one can get the 6 + — > KN pseudoscalar coupling from the generalized Goldberger-Treiman relation 

g A (e^KN)(M e + M N ) 
geKN = = 2.24 (79) 

where we use M e = 1530 MeV, M N = 940 MeV, F K = 1.2F n = 112 McV. Knowing the transition pseudoscalar 
constant one can evaluate the + width from the general expression for the ^ + hyperon decay [34] 

r e = 2 • ff ^IPl - M f m K = 4 . 44 MeV (80) 

87T Mq 



where |p| = \/(M| - - m 2 K ) 2 - 4Mf J m 2 K /2M e = 254 MeV is the kaon momentum in the decay (m K = 
495 MeV), and we have put the factor 2 to account for the equal-probability K + n and K°p decays. 



XIII. THEORETICAL UNCERTAINTIES 



Unfortunately, in baryon physics we deal with a truly strong interaction case, meaning that all dimcnsionlcss 
quantities are generally speaking of the order of unity. There is no really small algebraic parameter in sight that 
would allow some kind of perturbative expansion. We have argued in the Introduction that 1 /N c can be considered 
as a formal small parameter justifying the use of the Relativistic Mean Field Approximation. However, it is definitely 
not small enough when it comes to "kinematical" factors related to the rotational states of the mean-field baryons. 
Therefore, we treat the octet, decuplet and antidecuplet baryons as it should be at iV c = 3, rather than dealing 
with the large- iV c prototypes of those multiplets. We expect that the accuracy of this mixed logic is at the level of 
1/(2ttN c ) ~ 10%. 

Another source of the uncertainty is the present lack of knowledge of the exact low-energy effective action (2), in 
particular of the exact dynamical quark mass M{p). We have mimicked the fall-off of this function at large momenta 
by introducing the Pauli-Villars cutoff such that the F v constant and the chiral condensate <qq> are reproduced. 
From the experience in calculating various observables in the Chiral Quark Soliton Model [3] we estimate the ensuing 
error as <~ 15%. Thus, the resulting accuracy of the Relativistic Mean Field Approximation with exact account for 
the rotational wave functions of baryons, is expected to be about 20%, and this is indeed the typical accuracy with 
which formfactors, magnetic moments, parton distributions etc. have been computed in the model; in many cases the 
accuracy is actually much better but we quote here the pessimistically expected accuracy. 

When dealing with hyperons containing strange quarks, one has to decide how does one treat the mass m s . Theo- 
retically speaking, there are two small parameters, l/N c and m s /A where A is the characteristic scale of the strong 
interactions. Before choosing a calculational scheme one has to decide which of the two parameters is "smaller". One 
observes that the mass splittings in the baryon octet and decuplet are 0{m s N c ) and are somewhat less than the 
splitting between octet and decuplet centers, which is 0(A/N c ). Also, the Gell-Mann-Okubo relations arc satisfied to 
the 0.5% accuracy, which can be algebraically written as 0(m 2 /A 2 ). It indicates that the former parameter is larger 
than the latter, moreover it is not unreasonable to say that the strange quark mass is very small, m s — 0(A/N 2 ). 
In practical terms it means that in baryons, m s can be treated as a perturbation in most cases. In this paper, we 
have actually used the chiral limit, m s = 0, i.e. the zeroth order of that perturbation series. Computing first-order 
corrections in m s to the observables does not cause serious difficulties, see e.g. Refs. [28, 35], but we have not done 
it here. The penalty is expected at the 20% level. 

Within the Relativistic Mean Field Approximation, there arises a new important dimensionless parameter, namely 
e = E\ cv /M where M = M(0) is the dynamical quark mass at zero virtuality and -Ei cv is the quark bound-state level 
generated by the self-consistent chiral field, see eq. (6). If e w 1, the valence quarks in baryons are non-relativistic, the 
upper s-wave Dirac component of their wave function h(r) is much larger than the lower p-wave component j(r), and 
the number of additional QQ pairs in baryons tends to zero. In this limit the 6 + width goes to zero [35], which can be 
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also seen from the equations of the previous section, in particular from eq. (78): the numerator in that equation (A^) 
is proportional to the number of the QQ pairs while the denominator (y]VW(0)) is proportional to its square root. 
Consequently, the width Tq is proportional to the number of QQ pairs in ordinary baryons and vanishes in 
the non-relativistic limit e — > 1. 

Actually in our estimates in Section XII, we have systematically used the first-order perturbation theory in the 
"relativism" of valence quarks or, mathematically speaking, in 1 — e. Namely, we have 

ignored the lower component of the valence wave function j(r) 

ignored the distortion of the valence wave function by the sea, eq. (36) 

used the approximate expression for the pair wave function (28) 

computed the direct but neglected the exchange diagrams when evaluating the 5Q normalization and transition 
matrix elements, shown in Fig. 7 and 9 



• neglected the 7Q,9Q... components in baryons. 

It is difficult to evaluate the errors of these approximations before the neglected corrections are computed (which is 
surely feasible as all corrections are well defined above, but it has not been done). Unfortunately, the uncertainty 
associated with this non-relativistic approximation is expected to be large since the actual expansion parameter 
1 — e = 0.42 is poor. Another sign that the nuclcon is in fact a relativistic system comes from the 50% ratio of the 
5Q to the 3Q normalization. Treating the relativistic system in the first order in the "relativism" , is undoubtedly the 
main source of the uncertainty in our numerical estimates. 

Assuming that the uncertainties mentioned above are "statistically independent" , we estimate the error in computing 
the transition coupling gexN as 

-\/0.2 2 + 0.2 2 + 0.42 2 = 0.5 

implying a 100% error in the width. 

To get a feeling of the accuracy of our estimates, we have redone the calculations of Section XII replacing the 
probability distribution $(z,qjj introduced in Section XI by a flat one. This is a legitimate assumption within 
the Mean Field Approximation as it corresponds to ignoring the restriction following from the quark momentum 
conservation. We remind the reader that we have used the value of the baryon mass A4 = 1207 MeV instead of, say, 
940 MeV: the difference is believed to be partially due to adding the momentum conservation correction to the Mean 
Field result [36]. Therefore, it may seem to be more logical to ignore the quark momentum conservation systematically 
throughout the calculations. 

With this assumption, the evaluation of the 5Q matrix elements (59-62) is very easy and we obtain, instead of 
eq. (68), 

K nir = 0.0428, K aa = 0.0235, K 33 = 0.02 14, K 3lJ = 0.0226. (81) 
These numbers lead, via eqs. (64-67), to the following values of the physical quantities: 

-^^= 0.405, gf(N) = l.U, (82) 

which are not qualitatively different from the estimates (76,77). However, the 6 + width appears to be quite sensitive 
to the change: 

g A (& -> KN) = 0.146, T e = 2.32 MeV, (83) 

the width turning out nearly twice smaller than that of Section XII. It gives the idea of the accuracy of our estimate. 

Probably the worse error in our estimate of the + width arises from neglecting the exchange diagrams in matrix 
elements, see Figs. 7 and 9. As a rule, their account in processes involving fermions reduces matrix elements. It should 
be also noticed that the mass difference between the + and the nucleon is not small whereas we have estimated the 
transition amplitude at zero momentum transfer. One would hence expect that there is an additional formfactor-like 
reduction of the 9 + — ► KN transition amplitude. 

Therefore, one can well imagine that the 6 + width (83) is further reduced, maybe even below the 1 MeV value. We 
do not think that taking into account the 7Q... components in the transition matrix elements will seriously alter the 
5Q estimates. 

Pinning down the <d + width even inside a wide 50% error margin requires much more work than presented here. 
Nevertheless, the estimate that Tq is in a few MeV range seems to be safe. It follows from the relative suppression 
of QQ pairs in the nucleon, and from the SU (3) group suppression in the + — ► KN transition. 
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XIV. CONCLUSIONS 



Ordinary baryons are not made of three quarks only but have a substantial component with additional QQ pairs. 
For some observables, additional QQ pairs change the naive 3Q results by only 20% (like in the case of the nucleon 
axial constant) but for some other observables they change the naive result by a factor of 3 — 4 (as in the case of 
the spin carried by quarks or the nucleon a term). Hence it is imperative to learn how to work with higher Fock 
components in baryons. 

It is imperative not only for practical but for simple theoretical reasons. Assuming there are just three quarks 
in a baryon and wishing to write down their wave function, one realizes that one cannot "measure" (and hence 
mathematically describe) the quark position to an accuracy better than the Compton wave length of a pion (1.4 fm), 
since by uncertainty principle one then produces a pion or an additional QQ pair. Since the baryon size is 1 fm, there 
is literally no room for the just-three-quarks description of a baryon. The uncertainty principle demands that baryons 
should be described as containing an indefinite number of QQ pairs. The only question is quantitative: how many 
are there QQ pairs, and what are their wave functions [40]. 

Moving to this uncharted territory, one has to satisfy certain general conditions as the relativistic invariance (since 
pair production is a relativistic effect) and the completeness of states, needed to guarantee that parton distributions, 
including antiquarks, are positive-definite and are subject to sum rules following from the conservation laws for the 
baryon charge, axial current, etc. Relativistic invariance and the completeness of states can be achieved only in a 
relativistic quantum field theory. A field-theoretic model of baryons, which takes into account the infinite number of 
degrees of freedom and in which these general conditions are automatically met, is the Chiral Quark Soliton Model [3], 
an alias for the Relativistic Mean Field Approximation. 

Using this model, we have presented a technique allowing to write down explicitly the 3Q, 5Q, 7Q... wave functions 
of the octet, decuplet and antidccuplct baryons. In the exotic antidecuplet the 3Q component is, of course, absent, 
but its leading 5Q component is space-wise similar to the non-leading 5Q component of the nucleon. The technique is 
mathematically equivalent to the "valence quarks plus Dirac continuum" method exploited previously, but brings the 
mean field approach even closer to the language of the quark wave functions used by many people. We have shown that 
the standard SU (6) wave functions are easily reproduced for the 3Q components of the octet and decuplet baryons, if 
one assumes the non-relativistic limit. However, we have given explicit formulae for the relativistic corrections to the 
3Q wave function, and also for the 5Q wave function of the nucleon and of the exotic + . Having patience one can 
go further and write down e.g. the 19-quark component of the proton or the 7-quark component of the exotic S . 

It is important that the QQ pair in the 5Q Fock component of any baryon, be it the nucleon or the 9 + , is added 
in the form of a chiral field, which costs little energy. This is the reason why the 5Q component in the nucleon turns 
out to be substantial, and why the exotic 6 + baryon whose Fock decomposition starts from the 5Q component, is 
expected to be light. The energy penalty for making a pentaquark would be exactly zero in the chiral limit and 
were baryons infinitely large. In reality, to make e.g. the 6 + from the nucleon, one has to create a quasi-Goldstone 
K-meson and to confine it inside the baryon of the size > 1/M. It costs roughly 



Therefore, one should expect the exotic 9 + around 1540 MeV. The existence of the lightest degree of freedom in 
QCD, namely the pseudo-Goldstone fields, is ignored in the non-relativistic constituent quark models, which leads to 
the overestimate of the G + mass by typically 500 MeV [19]. 

Having presented the general formalism for computing observables for the 3Q as well as for higher Fock components, 
we have applied it to several cases of immediate interest. We have estimated the normalization of the 5Q component 
of the nucleon as about 50% of the 3Q component, meaning that about 1/3 of the time the nucleon is "made of five 
quarks. We have also shown that the account for the 5Q component in the nucleon moves its axial charge from the 
naive non-relativistic value of 5/3 much closer to the experimental value. 

Another case of interest is the width of the exotic + baryon: if it exists, why is it so narrow? The best direct 
experimental limit is Tq < 9 MeV [37], however indirect estimates [38] indicate that the width can be as small as 
1 Mcv or even less. Such a narrow width for a strongly decaying baryon some 100 MeV above the threshold, is the 
main surprise about the + . Since the original narrow-width estimate Tq < 15 MeV [35] (or, to be more precise, 
3.6 < Tq < 11.4 MeV [39]) based on the Chiral Quark Soliton Model, we have made here the first estimate of the 
axial constant for the + — ► KN transition, based on the direct computation of the 5Q matrix element within the 
same logic. We have shown that the + width is proportional to the number of QQ pairs in nucleons and is thus 
naturally suppressed, as compared to the expected widths of baryons with the dominant 3Q component. Assuming 
the SU(3) symmetry, the 9 + width is additionally suppressed by the SU(3) Clebsch-Gordan factors. 

In this first direct estimate using the 5Q wave functions of the + and of the nucleon, we have made several 
approximations summarized in Section XIII. The worse approximations can be eliminated by further work outlined 




(84) 
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in the paper but at the moment they lead to a large theoretical uncertainty in the + width. Depending on the 
way we impose the approximation, we obtain Tq w 2 — 4 MeV, with a high probability that it is further reduced by 
taking into account the quark exchange processes in the 9 + — > KN transition, and the formfactor-like suppression in 
this finite momentum transfer decay (both of which we neglected). Therefore, the 6 + width of a few MeV appears 
naturally within the Rclativistic Mean Field Approximation, without any parameter fixing. 

We believe that the presented formalism has a broad field of applications, apart from exotic baryons. One kind of 
applications has been already started in Ref. [16] and involves exclusive processes, nucleon distribution amplitudes, 
parton distributions for a fixed number of quarks, and the like. Another kind of applications is for low energies. 
One can compute any type of transition amplitudes between various Fock components of baryons, including the 
relativistic effects, the effects of the SU{3) symmetry violation, mixing of multiplcts, and so on. 
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APPENDIX A: PARAMETRIZATION OF SU(N) MATRICES 

In this Appendix we construct by induction a parametrization of a general unitary unimodular SU(N) matrix in 
terms of N 2 — 1 "Euler angles", and write down the invariant Haar measure over the group in terms of these angles. 
The construction has been prompted by the parametrization of the SU(3) group by Mathur and Sen [41]. 

The idea is to write iteratively a general SU(N) matrix as 

Rn = SnRn-i (Al) 

where Rn-i is a general SU (N—l) matrix with (N— l) 2 — 1 parameters and Sn is an SU(N) matrix of a special kind 
with only 2N — 1 parameters belonging to the sphere S' 2Ar_1 . It gives the full parametrization of a general SU(N) 
matrix with N 2 — 1 parameters. Accordingly, the invariant integration measure over the SU(N) group is presented 
as a product of measures over the spheres S 3 x S 5 x ... x S 2N ~ X . 

One starts from the SU{2) group whose parametrization as a 3d sphere S 3 is well known: one writes a general 
SU(2) matrix as 

a _( e~ mil cos0i e M12 sin0i\ , . 

2 ^ -e-^sin^! e Ml1 cos^ ) { > 

where the last column in S 2 can be viewed as a 2d complex vector v 2 = (z 1 , z 2 ) normalized as Iz 1 ! 2 + |z 2 | 2 = 1, which 
defines an S 3 sphere. The first column is the orthogonal vector v\ = e^vy. The group measure can be written as an 
integral over the S 3 sphere, 

-jjj J dz 1 dz 1 dz 2 dz 2 5(\z x \ 2 + \z 2 \ 2 - 1), (A3) 
or, explicitly in terms of three angles, as 

■ / d<fii sin <pi cos (j>\ / dan / dai 2 (= 1). (A4) 
Jo Jo Jo 



2tt 2 Jo 

To construct a general SU(3) matrix using the recipe (Al) we first make a 3 x 3 matrix i? 2 putting S2, say, in its 
left upper corner, 

R 2 =\ S2 ) , (A5) 
V 1 J 

and define 

e M23 cos 9 e M23 sin ( 

Rz = S 3 R 2 , S3 = I -e M22 sin 6 sin (j) 2 e~ iQ21 - iQ23 cos 2 e iQ22 cos 6 sin 2 |- (A6) 

. e »a 2 i sin#cos02 — e _M22_M23 sin </) 2 e M21 cos 8 cos 62 
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The last column in £3 can be viewed as a 3d complex vector V3 = (z 1 , z 2 , z 3 ) normalized to \z 1 \ 2 + \z 2 \ 2 + \z 3 \ 2 = 1, 
which defines an S 5 sphere. The three columns are constructed as (complexified) orts in spherical coordinates: 
vi <~ e r , v 2 ~ 60, i; 3 ~ eg. There is of course a freedom of choosing the orts and the angles; we use part of this freedom 
in such a way that R3 = I3 when all angles are set to zero. 
The measure on S 5 can be written as 

A J dz 1 dz 1 dz 2 dz 2 dz 3 dz 3 5(\z 1 \ 2 + \z 2 \ 2 + \z 3 \ 2 -l), (A7) 

or, explicitly in terms of five angles, as 

I r§ r§ r 2 ^ r 27T r 27T 

— / d9 cos Osm.6 I d(p 2 sin 2 cos 02 / da 2 i / da 2 2 / da 2 
n Jo Jo Jo Jo Jo 
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(= I)- (A8) 



The integrations limits are chosen such that the S 5 sphere is covered once. 

The full SU (3) measure is found in the standard way: one constructs the metric tensor 



dpm dp 

then the SU(3) measure is 



dR dR) 

Tr a — ~x~~^> P m = aii,oii 2 ,<j>i,a 2 i,a 22 ,a 2 3,<j) 2 ,9, m,n=1...8; (A9) 



\J det g ~ (sin 0i cos0i) • (cos 3 #sin#sin0 2 cos 02 ) (A10) 

i.e. it is factorized into the product of the measures over the spheres S 3 and S 5 , see eqs.(A4,A8). All group integrals 
in Appendix B can be performed directly using the above paramctrization of the SU(2) and SU(3) matrices and the 
above Haar measure. In fact we have checked the results of Appendix B in this way. 

The construction can be iteratively generalized to higher SU (N) groups in such a way that the group parameter 
space is a direct product of the spheres S 3 xS 5 x. . . S 2Ar_1 with the total number of parameters X^j=i (2J+1) = A 2 — 1, 
as it should be for the SU (N) group. For example, a parametrization of R £ SU (4) is 

Ra = SiR 3 , (All) 

( cosxe M34 -sinxe M34 \ 

sinxsin^e 1 " 33 cos^e 1033 cosxsin^e 1 " 33 

sin x cos 3 sin 03 e M32 - sin 9 3 sin 3 e M32 cos 3 e" M31 ~ M33 " M34 cosxcos6> 3 sin0 3 e~ M32 

y sinxcos6» 3 cos0 3 e M31 - sin 63 cos 3 e M31 - sin 3 e~ M32 ~ M33 " M34 cosxcos6» 3 cos0 3 e M31 J 

and i?3 is the block-diagonal 4x4 matrix with the general SU (3) matrix (A6) in the left upper corner and unity in the 
right lower corner. We thus add 7 new parameters to the previous 8 of the SU (3) parametrization. The columns of S4 
are complex orts in spherical coordinates, e r , eg, e$, e x . Denoting the last column v 4 = (z 1 , z 2 ,z 3 , z 4 ) the integration 
measure is that of a sphere S 7 : 

J dz 1 dz 1 . . . dz i dz i S^z 1 ] 2 + . . . |z 4 | 2 — 1) <~ J cos 5 xsinx cos 3 #3 sin #3 cos 03 sin0 3 dxd63d(f)3da3ida32da33da34. 

(A12) 

The full ST/ (4) measure built from the general rule (A9) (where now there are 15 angles) is factorized into the product 
of measures over the S 3 , S 5 and S 7 spheres. 

APPENDIX B: GROUP INTEGRALS 

In this Appendix, we give a list of group integrals used in the main text, over the Haar measure of the SU(N) 
group, normalized to unity, j dR — 1. 
For any SU (N) one has 

J dR R{ =0, j dR Ry = 0, J dR R{ R\i = Sf . (Bl) 
For N = 2 the following group integral is non-zero: 

J dRR{R^= l -e^e t3 . (B2) 
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For N > 2 this integral is zero; its analog in SU(3) is 

dRR[R]R h k ^ X -e^ h e l2k . (B3) 

On the contrary, in SU (2) it is zero. 

The general method of finding integrals of several matrices R, R* is as follows. The result of an integration over the 
invariant measure can be only invariant tensors which, for the SU (N) group, can be built solely from the Kronecker 
S and Levi-Civita e tensors. One constructs the supposed tensor of a given rank as a combination of S's and e's, 
satisfying the symmetry relations following from the integral in question. The indefinite coefficients in the combination 
are then found from contracting both sides with various <5's and e's and thus by reducing the integral to a previously 
derived one. 

For any SU (N) group one has 



dR R fl R^ jl R f2 R^ 32 — - S fl S f2 I (f. 1 ^ 2 —S h S i2 I +S fl S f2 I S h S i2 —S^S 1 - 2 

alttt h n gi n i 2 n 92 - N 2 _ i °gi°92 [°ji°j2 N°h°h J + °92°gi \°n°]i h°h 



(B4) 



since its contraction with, say, S 9 ^ must reduce it to eq. (Bl). 
In SU (2) there is an identity 

S i e hh + s i £j2h + s j 2 £j 3 ji =0, (B5) 
using which one finds that the following integral is non-zero: 

JdRR^Rf 2 Rf 3 R^ = \ (tf^eAAe^, + 6^ h ^e Ml + 6*S* h ^e hh ) . (B6) 

In SU(3) and higher groups this integral is zero. The analog of the identity (B5) in SU(3) is (notice the signs in the 
cyclic permutation!) 

^jl C 323334 ~ ^j 2 € 3334jl + ^j 3 e 343l32 ~ ^j 4 e 3l3233 = 0, (B7) 

and the analog of eq. (B6) is 

/ 

— — (s fl S 3 eMzft-e- ■ ■ +S f2 S j f^AAf ■ ■ ■ +S f3 (j 3 fhhhf ■ ■ + S f4 S 3 fhhhf- ■ ■ \ (BR) 

24 \ 9 31 t 323334 ^°g °j 2 t t 333i3i '°g u j 3 t 1 34.3132 ^ u g °j 4 t t 3i3233 J ■ \ DO ) 

This integral arises when one projects three quarks from the bound-state level onto the octet baryon. 
To evaluate the SU (3) average of six matrices, one needs the identities 

e hj233 e 3li2i3 i e i23233 e il3li3 e »3 32 33 e h *2 jl = 
e jlil33 e 32i2i3 i e 3li2j 3 e h32i3 ' e 3li 3 33 e il%233 = 

i e ili233 = e 3l3233 e ili2i3- (^9) 



dR R fl R{ 2 R{ 3 i?{ 4 Rl 3 

3l 32 33 34 g 



One gets 



\Jlj2U C J3*2*3 ~ c Ji;f2*2 c *lj3«3 ~ c Jlj2«3 c ili2j3 c Jl 32 33 c *l *2 *3 ' 



e iii2i3 e ii*2i3 



I f h 1 f 2 f3 e fih 2 h 3 ...... i JiihlaJiilihs ...... i ,'13/2/3 e hih 2 fi , . ...... 

T t t t »l.j2.j3 C .Jl*2»3 ' t E C »2j2j3 C »lJl»3 T t t t *3j2j3 t *l«2jl 

1 ,/i hif 3 f 2 h 2 h 3 ...... 1 f fih 2 f 3 e hif 2 h 3 . ...... 1 , /1 ^3/3 ,^1^2/2 , . ...... 

T C t t JlHj3 C j2»2J3 ~ fc E C Jl«2j3 C *1.J2»3 ~ t fc t Jl *3 J3 C H *2 J2 

i f fihhi foh2h 3 ...... i ,/1/2/M Jiif 3 h 3 ...... i f flf2h 3 hih 2 f 3 . . . . \ ('Rl 0^ 

T t t c Jl.j2»l t .j3»2l3 ~ fc t C JU2»2 C *1.J3*3 ~ t fc t Jl J2*3 E «l *2 J3 / ' V JJ±U ^ 

The result for the next integral is rather lengthy. We give it for the general SU(N). For abbreviation, we use the 
notation 

6 t 5 t 5 t 5 l S % S t = (231)(321), etc. (Bll) 
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One has 



/ 



dRRjl R-fj* Rjl -R^ 1 Rli 2 < 3 



N(N 2 - 1){N 2 -4) 

• {(TV 2 — 2) [(123)(123) + (132)(132) + (321)(321) + (213)(213) + (312)(231) + (231)(312)] 

- N [(123) ((132) + (321) + (213)) + (132) ((123) + (231) + (312)) + (321) ((312) + (123) + (231)) 
+ (213) ((231) + (312) + (123)) + (312) ((213) + (132) + (321)) + (231) ((321) + (213) + (132))] 
+ 2 [(123) ((312) + (231)) + (132) ((213) + (321)) + (321) ((132) + (213)) 

+ (213) ((321) + (132)) + (312) ((123) + (312)) + (231) ((231) + (123))]} . (B12) 

Apparently at N — 2 something gets wrong. For N — 2 there is a formal identity following from the fact that at 
N = 2 one has t ht * h €h lha h a = 0: 

(123) - (132) - (321) - (213) + (312) + (231) = 0. (B13) 
Consequently for SU(2) one obtains a shorter expression: 

dRR? 1 // ; // ; /<v : rV 2 rV 3 

31 32 33 hi tl2 tl3 

= - {[(123)(123) + (132)(132) + (321)(321) + (213)(213) + (312)(231) + (231)(312)1 
6 

- J [(123) ((132) + (321) + (213)) + (132) ((123) + (231) + (312)) + (321) ((312) + (123) + (231)) 
+ (213) ((231) + (312) + (123)) + (312) ((213) + (132) + (321)) + (231) ((321) + (213) + (132))]} . 

In case one is interested in the presence of an additional quark-antiquark pair in an octet baryon, one has to use 
the group integral 



360 





\ e e 3U2 [<>g Oft \ li0 j 4 °33 








+ 


Jlhh . 








^32) 


+ 


Jlhh . 

fc t 3l34 


H 2S fl K 5 ^ - 






^2) 


+ 


J2hh . 


- 








+ 


f hhh . 

fc fc j2j4 


K 5 ^ - 


^33) 






+ 


f hhh . 

C t ]334 











}• (B14) 

This tensor defines, in particular, the five-quark wave function of the nucleon, see eq. (45). 

For finding the quark structure of the antidecuplet, the following group integrals are relevant. The (conjugate) 
rotational wave function of the antidecuplet is (sec subsection IV. C) 



A*{hi h 2 h 3 } 1 ri\ _ 1 / n/ii r>h 2 r>h 3 , r/h 3 phi nh 2 , r/h 2 r>h 3 r>hi \ 
A k (K) - -(K 3 K 3 H k + K 3 H 3 K k + H 3 H 3 H k ). 



(B15) 



Projecting it on three quarks and using eq. (BIO) we get an identical zero because all terms in cq. (BIO) are antisym- 
metric in a pair of flavor indices while the tensor (B15) is symmetric. It reflects the fact that one cannot build an 
antidecuplet from three quarks: 



J dRRjl R f / 2 R f / 3 A* k {hl h2 ha} (R) = 0. 



(B16) 
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However, a similar group integral with an additional quark-antiquark pair is non-zero: 

J dR Rjl /,';;; //f (//;; Ft* ) T^ 1 h ° h *>(i?) = Jl_ {,,,,., ,, ; [^V 1 ' 2 '"^'" + ^h^hf^ 

+ tf 1 } Uhh^i e Sz}ih 3 _|_ e fif2h 3 e f 3 f4,h 2 \ _|_ fih 2 /f 1 f 2 h 1 ^f 3 fih 3 _|_ £ /i/ 2 e hfihi \ 



+ 



+ c,- 



§h 3 r e f2f3hi e fif4.h 2 _|_ e f 2 hh 2e hfihi\ _|_ ^ rf 2 f 3 h 2 e hfih 3 _|_ e f 2 f 3 h 3 e hfih 2 \ 



T e jlj 3 e 3234 



gh 3 ^ e f 1 f 3 h le f 2 f i h 2 _|_ e hf 3 h 2e hf4,h^ 
+ §^1 ^ e hhh 2e f 2 Uh 3 _|_ e fif 3 h 3e hfih 2 ^ _|_ fih 2 ^ e fif 3 h! e f 2 fih 3 _|_ e fif 3 h 3e hhhi | 

In particular, for the 6 + baryon being the 333-component of the antidccuplet we have 

6£0R) = V30^ 333 (i?) = V30i?li?li?^, 6 fc (i?) = V30i?f ^ 3 i?3, fe . 
The projection of five quarks onto to the + rotational wave function (B18) gives the tensor 



180 JU2 t j3j4 t fc I t j2j3 t JU4 t fc ' t JU3 t j2j4 t C ; 



(B17) 



(B18) 



(B19) 



This equation leads immediately to the five-quark wave function of the + , see eqs.(46,47). 
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